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I.  INTRODUCTION 

Consider  electromagnetic  coupling  through  a  small  aperture  in  an 
infinitely  thin  perfectly  conducting  screen  separating  two  regions.  Im¬ 
pressed  sources  in  one  or  both  regions  cause  electromagnetic  power  to  flow 
through  the  aperture.  The  field  due  to  the  aperture  is  the  actual  field 
minus  the  short-circuit  field.  The  short-circuit  field  is  the  field  that 
would  result  if  the  aperture  were  closed  with  an  infinitely  thin  sheet  of 
perfect  electric  conductor.  The  problem  is  to  find  the  far  field  due  to 
the  aperture. 

The  method  of  solution  is  similar  to  that  used  in  [1],  and  is  sum¬ 
marized  as  follows.  We  close  the  aperture  with  an  infinitely  thin  sheet 
of  perfect  electric  conductor  and  provide  for  the  tangential  electric 
field  originally  present  in  the  aperture  by  attaching  an  unknown  magnetic 
surface  current  M  to  one  side  of  the  sheet  and  -M  to  the  other  side  of  the 
sheet.  Continuity  of  the  tangential  magnetic  field  across  the  sheet  gives 
an  integral  equation  for  M.  This  integral  equation  is  then  solved  numeri¬ 
cally  by  means  of  the  method  of  moments.  The  solution  is  expressed  in 
terms  of  generalized  network  parameters  [2], 

Specifically,  M  is  expressed  as  a  linear  combination  of  the  quasi- 
static  magnetic  currents  M^,  and  M^,  which  are  equivalent  sources 
whose  far  fields  form  a  basis  for  the  far  field  of  M.  This  linear  combi¬ 
nation  is  substituted  into  the  integral  equation  for  M,  and  then  the 
integral  equation  for  M  is  tested  with  M^,  M^,  and  successively.  The 
result  is  a  system  of  three  equations  in  three  unknowns.  In  matrix  nota¬ 
tion,  these  equations  state  that  the  product  of  a  square  matrix  with  a 
column  vector  of  the  unknowns  is  equal  to  a  column  vector  of  known 
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coefficients.  This  matrix  equation  is  called  the  moment  equation.  The 
square  matrix  is  called  the  moment  matrix,  and  the  column  vector  of  known 
coefficients  is  called  the  excitation  vector. 

The  elements  of  the  excitation  vector  are  simply  related  to  the 
short-circuit  incident  fields  evaluated  at  the  aperture.  The  imaginary 
parts  of  the  elements  of  the  moment  matrix  are  obtained  by  equating  the 
solution  of  the  moment  matrix  to  the  solution  predicted  by  conventional 
small  aperture  theory  [3],  [4],  The  real  parts  of  the  elements  of  the 
moment  matrix  are  the  coefficients  appearing  in  the  quadratic  expression 
for  the  power  radiated  by  M. 

The  far  field  due  to  the  aperture  is  the  far  field  due  to  M.  The 
solution  to  the  moment  equation  gives  M  as  a  linear  combination  of  M  ,  My 
and  M^.  The  far  field  from  each  of  M^,  My  and  is  equal  to  the  far 
field  from  a  current  element. 

The  moment  solution  so  obtained  is  applied  to  three  examples.  The 
first  example  is  an  aperture  in  an  infinite  conducting  plane,  the  second 
is  an  aperture  in  a  transverse  wall  between  two  rectangular  waveguides, 
and  the  third  is  an  aperture  in  a  lateral  wall  of  a  rectangular  waveguide 
coupling  to  a  half-space.  In  all  three  examples,  the  far  field  (due  to 
the  aperture)  obtained  by  means  of  the  moment  solution  is  the  same  as  that 
obtained  by  means  of  the  radiation  reaction  method  devised  by  Collin  [5], [6] 

In  Appendix  A,  it  is  shown  that  the  moment  solution  satisfies  conser¬ 
vation  of  power.  Conservation  of  power  means  that  the  power  entering  the 
aperture  is  equal  to  the  power  leaving  the  aperture.  In  Appendix  B,  it  is 


shown  that  the  moment  solution  satisfies  reciprocity  in  the  sense  that 
remote  source  and  measurement  points  on  opposite  sides  of  the  aperture 
can  be  interchanged. 

II.  STATEMENT  OF  THE  PROBLEM 

Consider  a  small  aperture  in  an  infinitely  thin  perfectly  conducting 

plane  screen  separating  two  regions.  Region  a  to  the  left  of  the  aperture 

is  filled  with  loss-free  homogeneous  material  whose  permeability  is  y  and 

a 

whose  permittivity  is  £  .  Region  b  to  the  right  of  the  aperture  is  filled 

d 

with  a  loss-free  homogeneous  material  characterized  by  y,  and  e,  .  As  shown 

b  b 

ia  ia 

in  Fig.  1,  electric  and  magnetic  current  sources  J  and  M  are  impressed 
in  region  a,  and  sources  J[ib  and  Mib  are  impressed  in  region  b.  All  of 
these  sources  are  remote  from  the  aperture.  In  Fig.  1,  the  aperture  is 
labeled  A. 

The  problem  is  to  find  the  far  field  due  to  the  aperture  for  the 
situation  described  in  the  previous  paragraph.  By  definition  the  field 
due  to  the  aperture  is  the  actual  field  minus  the  short-circuit  field.  The 
short-circuit  field  is  the  field  that  would  result  if  the  aperture  were 
closed  with  a  plane  sheet  of  perfect  electric  conductor. 

In  Sections  III-VI,  a  solution  for  the  far-field  due  to  the  aperture 
is  constructed.  This  solution  is  valid  for  the  simple  situation  described 
in  the  second  from  the  last  paragraph.  However,  this  solution  can  serve  as 
an  approximate  solution  for  the  more  general  situation  in  which  the  screen 
is  not  plane,  inhomogeneities  such  as  conducting  walls  or  dielectric 
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obstacles  are  present,  and  power  Is  dissipated.  It  can  be  expected  that 
good  accuracy  will  be  obtained  if  the  screen  has  only  slight  curvature 
near  the  aperture,  no  inhomogeneities  occur  near  the  aperture,  and  no 
appreciable  power  is  dissipated  near  the  aperture. 

III.  DEVELOPMENT  OF  THE  MOMENT  SOLUTION 

The  electromagnetic  field  in  region  a  in  Fig.  1  is  simulated  by  the 

in  is 

situation  shown  in  Fig.  2.  In  Fig.  2,  the  original  sources  J[  and  M 
exist  in  region  a,  the  aperture  is  closed,  and  a  sheet  of  magnetic  current 
-M  is  placed  just  to  the  left  of  the  closed  aperture.  Here, 


M  =  £  *  n  (1) 

where  E  is  the  electric  field  in  the  aperture  in  Fig.  1  and  n  is  a  unit 

vector  normal  to  A.  As  shown  in  Fig.  2,  n  is  directed  away  from  region  a. 

In  Fig.  2,  (E*a,  jlia)  is  the  electromagnetic  field  due  to  (J^a,  M*a)  radi- 
ating  in  the  presence  of  the  complete  screen,  and  (E  (-M),  H  (-M))  is  that 
due  to  -M  radiating  in  the  presence  of  the  complete  screen.  The  complete 
screen  is  the  screen  with  the  aperture  closed.  The  total  field 
(Eia  +  Ea(-M),  Hia  +  Ha(-M))  in  region  a  in  Fig.  2  is  the  field  in  region 
a  in  Fig.  1. 

The  field  in  region  b  in  Fig.  1  is  simulated  by  the  situation  shown 
in  Fig.  3.  In  Fig.  3,  the  original  sources  J ^  and  Mib  exist  in  region  b, 

the  aperture  is  closed  ,  and  the  sheet  of  magnetic  current  M  is  placed  just 

to  the  right  of  the  closed  aperture.  In  Fig.  3,  (E*b,  H*b)  is  the  field 
due  to  (J*b,  Mib)  radiating  in  the  presence  of  the  complete  screen,  and 
(Eb(M),  Hb(M))  is  that  due  to  M  radiating  in  the  presence  of  the  complete 
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Fig.  2.  Equivalence  for  Region  a. 
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Elb  +  Eb(M),  Hlb  +  Hb(M) 


Fig.  3.  Equivalence  for  Region  b. 
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screen.  The  total  field  (E*b  +  _Eb  (M) ,  Hib  +  H^(M))  in  region  b  in  Fig.  3 
is  the  field  in  region  b  in  Fig.  1. 

The  sign  difference  between  the  magnetic  currents  in  Figs.  2  and  3 
is  dictated  by  continuity  of  the  tangential  electric  field  across  A  in 
Fig.  1.  Continuity  of  the  tangential  magnetic  field  across  A  in  Fig.  1 
requires  that 


-H*  (M)  -  (M)  =  H*b  -  H*3 

—tan  —  —tan  —  — tan  — tan 


in 


(2) 


where  the  subscript  tan  denotes  components  tangential  to  A.  Since  the 
aperture  is  small,  M  may  be  expressed  by  the  three  term  expansion 


M  =  +  V2M2  +  V£  (3) 

where  is  the  quasi-static  current  whose  radiation  field  is  that  of  a  unit 
magnetic  current  element  K£  ■  1  in  the  _t  direction,  M2  is  the  quasi-static 
current  whose  radiation  field  is  that  of  a  unit  magnetic  current  element 
K£  =  1  in  the  jt2  direction,  and  is  the  quasi-static  current  whose  radi¬ 
ation  field  in  region  a  is  that  of  an  electric  current  element  IS.  =  -juie^ 
in  the  n  direction.  The  t_  and  t_2  directions  are  orthogonal  directions 
in  the  plane  of  the  aperture.  All  of  the  above  three  current  elements  are 
placed  at  the  origin  of  coordinates.  This  origin  is  chosen  to  be  somewhere 
in  the  aperture. 

The  phase  of  each  of  M  ,  M2>  and  is  taken  to  be  constant  throughout 

the  aperture.  Because  their  radiation  fields  are  those  of  unit  magnetic 

current  elements,  and  M2  are  purely  real.  With  regard  to  M^,  the 

electric  current  element  IS.  =  -jwe  in  the  n  direction  is  equivalent  in 

a  - 
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region  a  to  an  Infinitesimal  loop  of  magnetic  current  In  the  plane  of 
the  aperture  characterized  by  [7,  p.  135] 

KS  =  1  (4) 

where  K  is  the  filamentary  magnetic  current  and  S  is  the  area  of  the 
loop.  It  follows  from  (4)  that  is  also  purely  real. 

To  obtain  the  moment  solution  to  (2),  we  substitute  (3)  into  (2), 
define  the  symmetric  product  <A,  IJ>  between  two  vectors  A  and  B  to  be 
the  integral  of  their  dot  product  over  A,  and  then  tak*-  the  symmetric 
product  of  (2)  with  M^,  •  and  M3  successively.  In  matrix  form,  the 

result  is 

[Ya  +  Yb]V  =  I1  (5) 

where  [Ya]  and  [Yb]  are  square  matrices,  and  V  and  I  are  column  vectors. 

3  b  i 

The  elements  of  [Y  ] ,  [Y  ] ,  and  I  are  given  by 


The  elements  of  V  are  the  coefficients  V^,  V^,  and  appearing  in  (3). 

The  elements  (6)  -  (8)  cannot  be  evaluated  in  a  straightforward 
manner  because  the  functional  forms  of  the  magnetic  currents  M^,  M^,  and 
are  not  known.  However,  it  has  been  established  that 
1)  The  magnetic  currents  M^,  M^,  and  are  real. 
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2)  The  radiation  field  of  is  that  of  the  unit  magnetic 
current  element  K2  =  1  in  the  _t  direction. 

3)  The  radiation  field  of  is  that  of  the  unit  magnetic 
current  element  K£  =  1  in  the  t_  direction. 

4)  In  region  a,  the  radiation  field  of  is  that  of  the 
electric  current  element  12,  =  -jwe  in  the  n  direction. 


The  magnetic  current  as  it  appears  in  H  (M^>  in  (6)  is  attached 
to  the  left-hand  side  of  the  complete  screen.  The  magnetic  current  M, 
as  it  appears  in  Hb(M^)  in  (7)  is  attached  to  the  right-hand  side  of  the 


complete  screen.  The  radiation  fields  of  and  are  those  of  the 
unit  magnetic  current  elements  mentioned  in  assumptions  2)  and  3)  re¬ 
gardless  of  the  side  of  the  complete  screen  to  which  and  M2  are 
attached.  Assumption  4)  states  that  when  is  attached  to  the  left- 

hand  side  of  the  screen,  its  radiation  field  is  that  of  the  electric 

current  element  12,  =  -jwe  in  the  n  direction.  However,  it  follows 

a  - 

from  the  result  of  [7,  prob.  3-6.]  that  when  is  attached  to  the  right- 
hand  side  of  the  complete  screen,  its  radiation  field  is  that  of  the 
electric  current  element  12,  =  -jwe,  in  the  n  direction. 

D  — 

IV.  EVALUATION  OF  THE  EXCITATION  VECTOR 


Assumptions  2)  to  4)  of  Section  III  will  now  be  used  to  evaluate 
the  elements  (8)  of  the  excitation  vector  I  .  In  view  of  the  definition 
of  the  symmetric  product,  (8)  becomes 


I 


i 

i 


(9) 


where 


10 


Hiads 


(10) 


rr 


A 


„ib. 

H  ds 


(ID 


Here,  ds  is  the  differential  element  of  area.  Application  of  reciprocity 
[7,  Sec.  3-8.]  to  the  right-hand  side  of  (10)  gives 


(-  jla  .  Ea(M.)  +  Mia  •  Ha(M.))ds 
—  —  — 1  —  — i 


(12) 


where  (Ea(M.),  Ha(M. ))  is  the  field  radiated  by  M,  attached  to  the  left- 
hand  side  of  the  complete  screen.  The  integration  in  (12)  is  over  the 

j.cl 

domain  of  2  and  M  .  Assuming  that  this  domain  is  in  the  radiation 
field  of  M^,  we  may,  by  virtue  of  assumptions  2)  to  4),  replace  (12)  by 


(-  Jia  •  E3 (M±)  +  Mia  •  Ha(M±))ds 


(13) 


where  is  the  unit  magnetic  current  element  KJ£  *  1  in  the  direction, 

$2  is  the  unit  magnetic  current  element  K2  =  1  in  the  direction,  and 

M,  is  the  electric  current  element  IS  =  -jioe  in  the  n  direction.  We 
—3  a  — 

emphasize  that  is  an  electric  current  rather  than  a  magnetic  current. 
The  impulsive  currents  M^,  and  M3  are  located  at  the  origin  of 
coordinates.  A9  it  appears  in  (13),  is  attached  to  the  left-hand 
side  of  the  complete  screen. 

Application  of  reciprocity  to  the  right-hand  side  of  (13)  gives 


ds  , 


i  =  1,2 


Eiads 


(14) 

(15) 


J 
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Because  of  the  impulsive  nature  of  M^,  M^,  and  M^,  (14)  and  (15)  reduce 
to 


Tla  -  Hla 

h  ~  Hi 


ia  _  ia 

h  "  H2 


(16) 


T3a  '  j“,E3a 

i3-  •]  ^  iti 

where  is  the  t^  component  of  at  the  origin,  is  the  t^  com- 

i-di  i.  cX  *| 

ponent  of  H  at  the  origin,  and  is  the  ri  component  of  E  at  the 
origin. 


A  development  similar  to  (12)  -  (16)  reduces  (11)  to 


ib  _  „ib 
1  "  H1 

ib  _  ib 
2  "  H2 


4 

j^b  e3 


(17) 


,  T,ib  ib  ,  „ib  .  ,  .  ,.ib 

where  ,  Hj  ,  and  are,  respectively,  the  _t  component  of  H  ,  the 

t^  component  of  and  the  n  component  of  E^3,  an  evaluated  at  the 

origin.  The  factor  in  the  last  of  equations  (17)  is  due  to  the  fact 

that  when  is  attached  to  the  right-hand  side  of  the  complete  screen, 

its  radiation  field  is  that  of  the  electric  current  element  Id  ■  -iwe.  in 

b 

the  n  direction  rather  than  Id  =  -jwe  in  the  11  direction.  Substitution 

Si 

of  (16)  and  (17)  into  (9)  gives 


I 

I 

I 


i 

1 

1 

2 

i 

3 


h?  - 


„lb  ia 
H2  ~  H2 


><SE3b  '  E.E3a> 


(18) 
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This  completes  evaluation  of  the  elements  (8)  of  the  excitation  vector  I  . 


V.  EVALUATION  OF  THE  IMAGINARY  PART  OF  THE  MOMENT  MATRIX 


3  b 

The  imaginary  parts  of  of  (6)  and  Y^  of  (7)  will  now  be 
evaluated  by  using  Be the' s  theory  [3].  According  to  the  extension  of  this 
theory  to  allow  for  contrasting  media  [4],  the  far  field  of  the  magnetic 
current  M  in  Fig.  2  can  be  approximated  by  the  far  field  of  the  combina¬ 
tion  of  the  magnetic  pole  dipole 


b  r  .  Tfib.  /..<Ld  ,,ibN  >• 

^m  =  y  +y  ["-l  aml(Hl  "  H1  *  "  -2  °m2(H2  "  H2 

a  b 


(19) 


and  the  electric  charge  dipole 


2e 


e  +e, 
a  b 


n  a 
—  e 


(e.E3* 


%E3b> 


(20) 


Furthermore,  the  far  field  of  the  magnetic  current  M  in  Fig.  3  can  be 

approximated  by  the  far  field  of  the  combination  of  the  magnetic  pole 

pa  £b 

dipole  — 2^  and  the  electric  charge  dipole  —  In  (19)  and  (20), 

^b  a 

a  ,  and  a  0  are  the  magnetic  polarizabilities,  and  a  is  the  electric 
mi.  nu  6 

polarizability.  The  polarizabilities  a  ,,  a  _,  and  a  are  purely  real 

ml  mz  e 

and  correspond,  respectively,  to  the  amyy*  an<*  ae  use<*  ln  f ® ] * 

Equation  (19)  is  simpler  than  [4,  Eq.  (34)]  which  contains  a  real 
symmetric  polarizability  dyadic.  However,  since  it  is  real  and  symmetric, 
this  dyadic  can  be  diagonalized  by  a  simple  rotation  the  coordinate  axes 
in  the  plane  of  the  aperture.  Hence,  (19),  in  which  the  magnetic  polari¬ 
zability  dyadic  is  represented  only  by  the  diagonal  elements  a  .  and  a  _, 

ml  nu 

can  be  obtained  for  any  small  aperture. 
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For  our  purposes,  It  is  more  convenient  to  deal  with  current 
elements,  magnetic  denoted  by  K2,  and  electric  denoted  by  1^.  These  are 
related  to  and  by 

KA  =  jwya  2^  (21) 


Ift  =  .jco 


(22) 


Thanks  to  (19)  and  (20),  (21)  and  (22)  become 
2jioy  y,  . 

—  =  7-TiT  ^  II  “mX*  -  H1 b>  -  «m2<H2  -  «2  >1  <23> 

a  b 

2jwe  .  ,. 

u  -  r-j3-  n  ae  <£  e“  -  )  (24) 

a  b 

The  combination  of  -Kit  and  -lit  replaces  -M  in  Fig.  2. 

Previously,  the  magnetic  current  M  in  Fig.  3  was  likened  to  the 

combination  of  the  magnetic  pole  dipole  —  and  the  electric  charge 
%  b 

dipole  —  £e<  These  dipoles  are  equivalent  to  magnetic  and  electric  cur- 
a 

rent  elements  KJt'  and  lit*  given  by 


SI  *  t  £m> 

D 


iil  =  jw  (^  £e) 

a 

Comparison  of  (25)  with  (21)  gives 


(25) 


(26) 


Kit’  =  Kit 

Comparison  of  (26)  with  (22)  gives 

e. 

n'  =  -f-  n 

a 


(27) 


(28) 
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The  combination  of  K£*  and  !£'  replaces  M  in  Fig.  3. 

Replacement  of  -M  in  Fig.  2  by  the  combination  of  -lOt  and  -l£  and 
replacement  of  M  in  Fig.  3  by  the  combination  of  K£*  and  I£*  gives  the 
situation  shown  in  Fig.  4.  In  Fig.  4,  K£  is  the  magnetic  current  element 
(23),  and  _I£  is  the  electric  current  element  (24).  If  the  aperture  is 
sufficiently  small,  the  electromagnetic  field  in  Fig.  4  is  a  good  approxi¬ 
mation  to  that  in  Fig.  1  at  points  sufficiently  remote  from  the  aperture. 
Note  that  if  the  result  of  [7,  Prob.  3-6.]  is  used  to  replace  the  electric 

£b 

current  elements  -I£  and  — L£  in  Fig.  4  by  infinitesimal  loops  of  magnetic 

a 

current,  then  the  magnetic  current  attached  to  the  left-hand  side  of  the 
screen  in  Fig.  4  becomes  the  negative  of  the  magnetic  current  attached  to 
the  right-hand  side  of  the  screen  in  Fig.  4.  This  is  in  agreement  with  the 
fact  that  the  magnetic  current  attached  to  the  left-hand  side  of  the  screen 
in  Fig.  2  is  the  negative  of  the  magnetic  current  attached  to  the  right-hand 
side  of  the  screen  in  Fig.  3. 

The  combination  of  current  elements  -K£  and  -l£  in  Fig.  4  replaces 
-M  in  Fig.  2.  The  radiation  field  of  the  combination  of  K£  and  _l£  attached 
to  the  left-hand  side  of  the  complete  screen  is  (approximately)  equal  to  the 
radiation  field  of  M  attached  to  the  left-hand  side  of  the  complete  screen. 
The  combination  of  K£  and  _I£  is  called  M  and  is  expressed  as 

M  =  V, M.  +  V.M,  +  VH  (29) 

—  1—1  2.  Z  3 

where  and  E3  are  defined  just  after  (13).  Thanks  to  (23)  and  (24), 

the  V’s  are  given  by 


ria 


V 


ia 


Region  a 

Ua»  e 
a  a 


-I£ 


-mi 


r 


M; 


rib 


Region  b 

V  S 


U 


mi 


<L 


Complete  screen 


Fig.  4.  Situation  obtained  by  replacing  -M  in  Fig.  2  and 
M  in  Fig.  3  by  combinations  of  magnetic  and  elec¬ 
tric  current  elements. 
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0 

0 

2a 

_ e _ 

Jw(ea  +  £b) 

where  I*,  I**  and  I*  are  given  by  (18). 

Equations  (29)  and  (30)  were  oKta  .ed  by  using  Bethe's  theory  extended 
to  allow  for  contrasting  media.  This  theory  is  called  conventional  small 
aperture  theory.  For  calculation  of  the  far  field,  the  moment  solution  (3) 
is  equivalent  to  (29)  where  thst  t  are  the  elements  of  the  column  vector  V 
which  satisfies  the  moment  equation  (5).  Hence,  the  moment  solution  will  be 
equal  to  the  result  of  conventional  small  aperture  theory  if  the  elements  of 
the  solution  V  to  (5)  are  given  by  (30) .  By  forcing  the  elements  of  the  solu¬ 
tion  V  to  (5)  to  be  equal  to  (30)  for  arbitrary  values  of  I*,  I*,  and 
I*,  we  obtain  that  [Ya  +  Yb]  1  is  equal  to  the  square  matrix  in  (30).  This 
result  implies  that 

ya  *  % 

0 

0 

3l  b 

In  order  to  isolate  [Y  ]  from  [Y  ] ,  we  rewrite  (31)  as 


(32) 


where  k  is  the  wave  number  in  region  a,  n  is  the  intrinsic  impedance 

d  ci 

is  the  intrinsic 

impedance  in  region  b.  From  (6)  and  (7)  and  the  nature  of  the  magnetic 

^  cl 

field  operators  H  and  H  ,  we  expect  the  elements  of  [Y  ]  to  be  pro¬ 
portional  to  1/r)  and  the  elements  of  [Yb]  to  be  proportional  to  1/rj  . 
It  is  now  evident  from  (32)  that 


in  region  a,  is  the  wave  number  in  region  b,  and  ri 


[Ya] 


2jojy  a  . 

J  a  ml 


[Yb] 


2J^baml 


0 

1 


2^aam2 

0 


0 

1 


2>Vm2 


0 

0 

2ae 

0 

0 

jwe. 


2a 


e  J 


(33) 


(34) 


The  admittance  matrices  (33)  and  (34)  came  out  purely  imaginary 

because  no  radiation  terms  were  included  in  the  magnetic  pole  dipole  (19) 

3  b 

and  the  electric  charge  dipole  (20).  Now,  [Y  ]  and  [Y  ]  have  real  parts 
which,  although  usually  small  compared  to  the  imaginary  parts,  become 
important  if  the  imaginary  parts  are  cancelled  by  means  of  electromagnetic 


interaction. 
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VI.  EVALUATION  OF  THE  REAL  PART  OF  THE  MOMENT  MATRIX 

ci  b 

In  this  section,  the  real  parts  of  Y  and  Y  are  evaluated  and 
the  far  field  due  to  the  aperture  is  obtained  from  the  moment  solution 


for  M. 


cl  & 

According  to  (6),  the  real  part  of  Y^ ^ ,  denoted  by  Re(Y^),  is 


given  by 


Re(Ya  )  =  -  Re<M.,  Ha(M  )> 
ij  — i  —  ~j 


Unfortunately,  (35)  is  not  amenable  to  calculation  because  the  quasi¬ 
static  magnetic  currents  M.  and  M.  are  not  known. 

~i 

With  a  view  toward  obtaining  a  workable  formula  for  Re(Y^),  we 
consider  the  power  P  supplied  by  the  magnetic  current  M  of  (3)  attached 
to  the  left-hand  side  of  the  complete  screen.  If  M  is  the  only  source  in 
region  a,  then  P3  is  given  by  (7,  Eqs.  (1-63)  and  (1-64)] 

Pa  =  -  Re  <M*,  Ha(M)>  (36 

where  the  asterisk  denotes  complex  conjugate.  Because  M^,  fL,,  and 
are  real,  substitution  of  (3)  into  (36)  gives 

Pa  -  Re(V*  Y*V)  (37 


where  V  is  the  column  vector  of  the  coefficients  V^,  V^,  and  V^  in  (3). 

~  -►  3 

Also,  V  is  the  transpose  of  V,  and  Y  is  the  square  matrix  whose  elements 
are  defined  by  (6).  Equation  (37)  can  be  rewritten  as 

Pa  =  |  V*  [Y3  +  [Ya]*]V  (36 

From  reciprocity,  Y  is  a  symmetric  matrix  so  that  (38)  reduces  to 


Pa  -  V*[Re  Ya]V 
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g 

If  the  medium  in  region  a  is  loss- free,  then  P  is  also  equal 
to  the  integral  of  the  real  part  of  the  outward  component  of  the  complex 
Poynting  vector  over  a  surface  S  such  that  the  combination  of  S  and  the 
complete  screen  encloses  the  magnetic  current  M.  Hence, 

Pa  =  Re  j  Ea (M)  x  Ha* (M)  •  ds  (40) 

S 

where  the  magnitude  of  djs  is  the  differential  element  of  area  and  the 
direction  of  ds  is  that  of  the  normal  vector  which  points  outward  from  S. 

If  S  is  sufficiently  remote  from  M,  then  (40)  becomes 

Pa  =  Re  ([  Ea$)  x  Ha*(fl)  •  ds  (41) 

JJ  -  -  - 

S 

where  M  is  given  by  (29).  The  right-hand  side  of  (41)  is  the  power  radiated 
by  the  source  M  attached  to  the  left-hand  side  of  the  complete  screen.  Sub¬ 
stitution  of  (29)  into  (41)  gives 

Pa  =  v*[Ya]V  (42) 

where  Ya  is  a  square  matrix  whose  elements  are  given  by 

ff.  -  i  ff  [Ea(Mj  x  HS*  (M.  )  +  E3*  (M.  )  x  H®  (M.  )  ]  •  ds  (43) 

ij  2  jj  -  -  — i  -  ~i  — 

S 

The  magnetic  current  M  was  introduced  in  (36)  as  an  independent 
source.  Hence,  V  is  arbitrary  in  (39)  and  (42).  The  fact  that  (39)  is 
equal  to  (42)  for  all  conceivable  column  vectors  V  implies  that 

Re[Ya]  -  [Ya]  (44) 

According  to  (44),  knowledge  of  the  power  radiated  by  the  source  M  attached 
to  the  left-hand  side  of  the  complete  screen  determines  Re[Ya] .  Retracing 
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the  development  (36)  -  (44)  with  M  attached  to  the  right-hand  side  of 
the  complete  screen,  we  obtain 

Re[Yb]  =  [Yb] 

where  Y  is  a  square  matrix  whose  elements  are  given  by 


f  =  ± 
ij  2 


ff 


j  [ E° (fij )  x  Hb*(M^)  +  (M^)  x  HD"(M,)]  •  ds 


b*. 


,b* , 


(45) 


(46) 


According  to  (45),  knowledge  of  the  power  radiated  by  the  source  M  attached 
to  the  right-hand  side  of  the  complete  screen  determines  Re[Yb], 

Equations  (44)  and  (45)  call  for  replacement  of  (33)  and  (34)  by 


and 


[Ya] 


[Yb]  = 


1 

.  va 

va 

va 

2j(4iaaml 

+  Y11 

Y12 

Y  13 

Y3 

*21 

+  Y22 

Ya 

23 

Y3 

31 

Y3 

32 

jcoe 

a 

2a 

e 

1 

+  Yb 

Yb 

vb 

2jWVml 

+  Y11 

Y12 

Y13 

Yb 

1  t  Yb 

Yb 

Y21 

2^%am2  22 

Y23 

Yb 

Yb 

j^ 

Y31 

Y32 

2a  + 

33  J 


33 


(47) 


(48) 


As  stated  just  after  (38),  [Y  ]  is  a  symmetric  matrix.  Furthermore, 
[Yb]  is  also  a  symmetric  matrix.  For  the  situation  in  Fig.  4  where 
homogeneous  media  exists  on  both  sides  of  the  complete  screen,  [Ya] 
and  [Yb]  are  diagonal  matrices.  However,  the  elements  of  [Ya]  and  [Yb] 


as  obtained  from  (A3)  and  (46)  are  valid  for  the  more  general  situation 

in  which  the  complete  screen  is  not  plane  and  loss-free  inhomogeneities 

such  as  perfectly  conducting  walls  or  loss-free  dielectric  objects  are 

present  on  either  or  both  sides  of  the  complete  screen.  Presumably,  a 

slight  curvature  of  the  complete  screen  or  loss-free  inhomogeneities 

remote  from  the  aperture  will  have  little  effect  on  the  imaginary  parts 
a  b 

of  [Y  ]  and  [Y  ].  Hence,  it  is  feasible  to  apply  (47)  and  (48)  when  the 
complete  screen  curves  slightly  and/or  loss-free  inhomogeneities  exist 
remote  from  the  aperture.  In  these  cases,  [^a]  and  could  have 

non-zero  off-diagonal  elements. 

/-'ci  ~b 

The  matrix  elements  Y..  and  Y. .  in  ■‘'uressions  (47)  and  (48)  for 

ij  ij 

ci  b 

the  moment  matrices  [Y  ]  and  [Y  ]  are  defined  by  (43)  and  (46).  In 
(43),  is  the  unit  magnetic  current  element  Kf  =  1  in  the  t_^  direc¬ 
tion,  is  the  unit  magnetic  current  element  K(  3  1  in  the  t^  direc¬ 

tion,  and  is  the  electric  current  element  If.  =  -jtoe  in  the  n 
—3  a  — 

direction.  As  they  appear  in  (43),  M^,  M^,  and  are  attached  to  the 

left-hand  side  of  the  complete  screen.  In  (46),  Is  the  unit  magnetic 

current  element  Kf  =  1  in  the  direction,  is  the  unit  magnetic  cur¬ 
rent  element  Kf  =  1  in  the  direction,  and  is  the  electric  current 
element  If  in  the  n  direction.  As  they  appear  in  (46),  M^, 

and  are  attached  to  the  right-hand  side  of  the  complete  screen. 

The  magnetic  current  elements  and  in  (46)  are  the  same  as  those 
in  (43).  However,  the  electric  current  element  in  (46)  is  different 

from  that  in  (43)  if  e  j  £,  .  The  use  of  in  (43)  and  (46)  is  based 

a  b  — j 

on  the  global  definition  that  is  the  electric  current  element  which  is 
equivalent  to  the  infinitesimal  loop  of  magnetic  current  in  the  plane  of 
the  aperture  characterized  by  (4). 
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Afl  A^J 

Expressions  (43)  and  (46)  for  the  matrix  elements  and 
appearing  in  (47)  and  (48)  were  derived  for  loss-free  media  in  regions  a 
and  b.  However,  (43)  is  valid  so  long  as  there  is  no  loss  inside  S. 
Hence,  it  is  feasible  to  apply  (43)  if  there  is  some  loss  provided  that 
the  loss  in  the  vicinity  of  the  aperture  is  negligible.  Similarly,  (46) 
still  applies  if,  in  region  b,  loss  occurs  remote  from  the  aperture. 

as 

Expression  (43)  for  Y_  was  extracted  from  the  power  radiated 
by  the  source  M  attached  to  the  left-hand  side  of  the  complete  screen. 
Since  the  medium  inside  S  is  loss-free,  this  power  must  be  equal  to  the 

/N  A.  £ 

power  supplied  by  M.  An  alternative  expression  for  will  now  be 

A  A  -A 

extracted  from  the  power  supplied  by  M.  Since  and  are  purely  real 

/v 

magnetic  currents  and  is  a  purely  imaginary  electric  current,  the 
power  supplied  by  M  is  given  by 


P  = 


Re  <V  +  V2M2  ,  Ha(M)>  +  Re  <V^3  ,  Ea(M)> 


i49) 


Substitution  of  (29)  into  (49)  gives 


q  a-* 

P  =  Re(V  C  V) 


(50) 


where  Ca  is  the  square  matrix  whose  elements  are  given  by 


Cij  =  "  ^i’  , 


C^j  =  43.  Ea(Mj)> 


|i  =  1,2 

jj  -  1,2,3 

j  -  1,2,3 


(51) 


(52) 


Equation  (50)  can  be  rewritten  as 


P  =  V  Y  V 


(53) 


y3  =  |  ica  +  [can 


where 


(54) 
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The  matrix  Ya  defined  by  (54)  is  equal  to  the  matrix  Ya  appearing  in 

(42)  because  the  quadratic  forms  (42)  and  (53)  are  equal  to  each  other 

— ► 

for  arbitrary  V. 

g 

From  reciprocity,  C  is  a  symmetric  matrix  so  that  (54)  reduces  to 

Ya  =  Re  Ca  (55) 


Thanks  to  the  definition  of  the  symmetric  product  and  the  definitions 

A  A  A  ^ 

of  M^,  M,,,  and  M^,  substitution  of  (51)  and  (52)  for  the  elements  of  C 


in  (55)  gives 


~a 

Y  = 


-Re  Ha  (M^ )  -Re  Ha  (M 2  ) 

-Re  Hj^)  -Re  Ha(M2) 


cue  Im  Ea  (M  )  we  Im  Ea(M.) 
a  3  — 1  a  3  — 2 


-Re  H^(M3) 
-Re  Ha(Mj 


2  —3 
a 

3-3 


we  Im  Ea(M0)| 
a 


(56) 


where  Im  denotes  imaginary  part.  In  (56),  H^fM^),  ^(M^),  and  E^fftj)  are, 

g  A  g  A 

respectively,  the  _t^  component  of  H  (M^),  the  component  of  H  (M^),  and 

a  a 

the  n  component  of  E  (M^).  These  three  fields  are  evaluated  at  the 
origin.  In  (56),  M  is  the  electric  current  element  If,  =  -jwe  in  the  ri 

J  3 

direction. 

Consideration  of  the  power  supplied  by  M  attached  to  the  right-hand 
side  of  the  complete  screen  leads  to 


-Re  Hb(M:) 

-Re  Hb(M2) 

-Re  Hb(M^) 

Yb  = 

-Re  Hb (Mx ) 

-Re  Hb(M2) 

-Re  Hb(M3) 

(57) 

0)Cb Im  Eb(ML) 

■jc  Im  Eb(M2) 

i«Im  Eb(M3) 
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Here,  Hb(M^),  Hb(M..),  and  Eb(M.,)  are,  respectively,  the  component  of 

b  ])  a  b^ 

H  (M j ) ,  the  component  of  H  (M^),  and  the  n  component  of  E  (M^).  These 

three  fields  are  evaluated  at  the  origin.  In  (57),  is  the  electric 

current  element  l£  =  -jtoe,  in  the  n  direction. 

D  ‘ 

Now  that  I*,  Ya,  and  Yb  have  been  evaluated,  the  moment  equation 
(5)  can  be  solved  for  V.  The  solution  to  (5)  is 

-*■  a  h  —  1-M 

V  =  [Y  +  Y  ]  I  (58) 

~ a  b 

where  the  elements  of  I  are  given  by  (18).  The  elements  of  Y  and  Y 

are  given  by  (47)  and  (48)  in  which  Ya^  is  given  by  either  (43)  or  (56), 

^b 

and  Y  is  given  by  either  (46)  or  (57).  The  elements  V^,  V^,  and 
— ► 

of  V  of  (58)  give  M  according  to  (3). 

The  field  due  to  the  aperture  consists  of  the  fields  radiated  by 

-M  in  Fig.  2  and  M  in  Fig.  3.  If  only  the  far  field  is  of  interest,  M 

may  be  replaced  by  M  of  (29).  Hence,  in  region  a,  the  far  electric  field 

a  a 

due  to  the  aperture  is  _E  (-M)  given  by 

3 

Ea(-M)  =  -  l  V.  Ea(M  )  (59) 

i=l  1 

where  the  V's  are  the  elements  of  V  of  (58),  and  Ea(M. )  is  the  electric 

—  - 1 

A 

field  due  to  attached  to  the  left-hand  side  of  the  complete  screen. 

Here,  is  the  unit  magnetic  current  element  K£  *  1  in  the  _t  direction, 

is  the  unit  magnetic  current  k£  =  1  in  the  t^  direction,  and  is  the 

electric  current  element  I £  *=  -ju>e  in  the  n  direction. 

a 

Similarly,  in  region  b,  the  far  electric  field  due  to  the  aperture 
is  Eb(M)  given  by 

Ib(M)  -  l  V  Eb(M  ) 
i-1 


(60) 


— ►  Jj  A 

where  the  V's  are  the  elements  of  V  of  (58),  and  E  (M)  is  the  electric 

— i 

field  due  to  attached  to  the  right-hand  side  of  the  complete  screen. 
Here,  is  the  unit  magnetic  current  element  K£  =  1  in  the  t^  direction, 

A 

is  the  unit  magnetic  current  element  K£  =  1  in  the  t^  direction,  and 

M_  is  the  electric  current  element  I£  =  -jwe,  in  the  n  direction. 

— J  d  — 

VII.  APERTURE  IN  AN  INFINITE  CONDUCTING  PLANE 

In  this  section,  the  moment  solution  for  small  aperture  coupling 
is  applied  to  a  small  aperture  in  an  infinite  conducting  plane.  Here, 
the  situation  is  precisely  the  situation  described  in  the  first  para¬ 
graph  of  Section  II. 

In  expression  (58)  for  V,  the  elements  of  I*  are  given  by  (18), 

3  b 

Y  is  given  by  (47),  and  Y  is  given  by  (48).  The  matrix  elements  Y 
appear  in  (47),  and  the  matrix  elements  Y^  appear  in  (48). 

Now,  Y^j  appears  in  the  quadratic  expression  (42)  for  the  power 
radiated  by  the  source  M  of  (29)  attached  to  the  left-hand  side  of  the 
complete  conducting  plane.  Here,  the  adjective  complete  means  that  the 

A  A 

aperture  has  been  closed  with  a  perfect  conductor.  In  (29),  M^,  HL,, 
and  are  the  current  elements  defined  in  the  second  from  the  last 
paragraph  of  Section  VI.  In  the  case  of  the  complete  conducting  plane, 
there  are  no  cross-product  terras  in  the  quadratic  expression  for  the  power 

g  A  g 

P  radiated  by  M.  It  is  evident  that  P  is  twice  the  power  that  would 

A 

be  radiated  by  M  in  the  homogeneous  medium  characterized  by  (u  ,  E  ) 
everywhere.  From  [7,  Eq.  2-116]  for  the  power  radiated  by  an  electric 
current  element  in  a  homogeneous  medium  and  from  the  dual  expression  for 
the  power  radiated  by  a  magnetic  current  element  in  a  homogeneous  medium. 


we  obtain 


f 
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In  region  a,  the  far  electric  field  due  to  the  aperture  is  given 
by  (59)  in  which  the  V's  are  given  by  (64)  and  £a(M^)  is  the  far  electric 
field  due  to  attached  to  the  left-hand  side  of  the  complete  conducting 

A 

plane.  Here,  M^,  i=l,2,3,  are  the  current  elements  defined  in  the  second 
from  the  last  paragraph  in  Section  VI. 

It  is  apparent  from  [7,  Sec.  3-13]  that  the  far  electric  field 

£ 

E:  (If.)  due  to  an  electric  current  element  _If  attached  to  the  left-hand 
side  of  the  complete  conducting  plane  is  given  by 


a  Jt4Jae 

Ea(IJL)  =  -  a 


2iTr 


If. 


tan 


(65) 


where  r  is  the  distance  from  the  origin  and  If.  ^  is  the  part  of  Lf  per¬ 
pendicular  to  the  radius  vector  from  the  origin.  Similarly,  the  far 
magnetic  field  H  (Kf)  due  to  a  magnetic  current  element  Kf.  attached  to  the 
left-hand  side  of  the  complete  conducting  plane  is  given  by 


H  (Kf)  =  - 


-jk  r 
jtoe  e  a 


2?Tr 


Kf. 


tan 


(66) 


where  Kf.  is  the  part  of  Kf  perpendicular  to  the  radius  vector  from  the 
— tan  — 

origin.  In  the  far  zone,  the  electric  and  magnetic  fields  Ji  and  Jrf  due 
to  any  source  in  the  vicinity  of  the  origin  satisfy 


Ea  =  n  Ha  x  u 


(67) 


where  u  is  the  unit  vector  in  the  direction  of  the  radius  vector  from 
— r 

the  origin. 

Expressions  (65)  -  (67)  will  now  be  used  to  obtain  E!a(M^),  i  *  1,2,3. 
For  convenience,  an  xyz  coordinate  system  is  constructed  such  that  the  unit 
vectors  in  the  x,y,  and  z  directions  are  t^,  t^*  and  n,  respectively.  With 


this  coordinate  system,  (65)  gives 
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k2  e'jk*r  sin  6 
ia(M3)  -  Uq  - 


(68) 


where  0  is  the  angle  measured  from  the  positive  z  axis  and  is  the  unit 
vector  in  the  0  direction.  Similarly,  (66)  gives 

~jk  r 


jk  e 


H  (M3)  =  -  (u^  cos  0  cos  <j>  -  u^  sin  <}>)  ^  r 


and 


-jk  r 


jk  e 


£aQl2)  =  “  (Uq  cos  e  sin  $  +  cos  <t>)  ~ ~2n'r~T~ 


(69) 


(70) 


where  <t>  is  the  angle  measured  in  the  xy  plane  from  the  positive  x  axis 
and  u ,  is  the  unit  vector  in  the  <j>  direction.  Application  of  (67)  to 
(69)  and  (70)  gives 

-jk  r 
, _  J  a 

jkae 

£  (M^)  =  (jJq  sin  <t>  +  u^  cos  6  cos  <J>)  - 7^ - 


and 


•jk  r 


jk  e 

Ea(M„)  =  (-  u.  cos  <J>  +  u,  cos  0  sin  <J>)  3 

—  —2  —y  — $ 


2nr 


(71) 


(72) 


The  far  field  E°(M, )  due  to  M.  attached  to  the  right-hand  side  of  the 
—  — i  —1 

complete  conducting  plane  is  given  for  i  =  1,2,  and  3  by  equations  similar 
to  (71),  (72),  and  (68). 

..  -jv 

b  -  jkbe 

E  (M., )  =  (in,  sin  0  +  u,  cos  0  cos  <j>)  - -  (73) 

—  —1  “6 


2Ttr 


b  -  jkb" 

E  (M2)  =  (-  Uq  cos  p  +  u^  cos  0  sin  <J>)  —  27Tr 


*jV 


(74) 


l  k  e 

Eb(M.)  =  u 


2 


sin  0 


—3  -0  2Trr 

The  far  electric  field  Ea(-M)  in  region  a  due  to  the  aperture  is 


(75) 


given  by  (59)  which  is 


29 


3 

Ea(-M)  =  -  l  V  E3^)  (76) 

i=l 

where  the  V's  are  given  by  (64)  and  £a(M^),  i  =  1,2,3,  is  given  by  (71), 

b  A 

(72),  and  (68).  The  far  electric  field  E  (M)  in  region  b  due  to  the 
aperture  is  given  by  (60)  which  is 

Eb(M)  =  I  V.  Eb(M.)  (77) 

-  -  .  ,  1  -  - 1 

1=1 

where  the  V's  are  given  by  (64)  and  ^(M^),  i=l,2,3,  is  given  by  (73), 

(74),  and  (75). 

VIII.  APERTURE  IN  A  TRANSVERSE  WALL  BETWEEN  TWO  WAVEGUIDES 

In  this  section,  the  moment  solution  is  applied  to  two  rectangular 
waveguides  of  dimensions  a  by  b  and  a'  by  b'  separated  by  a  transverse 
wall  at  z  =  0  which  contains  a  small  aperture  as  shown  in  Fig.  5.  The 
aperture  is  located  at  the  center  (x  =  a/2,  y  =  b/2)  of  the  cross  section 
of  the  left-hand  waveguide  or,  equivalently,  at  the  center  (x'  =  a'/2, 
y'  =  b'/2)  of  the  cross  section  of  the  right-hand  waveguide.  Here,  x 
and  y  are  rectangular  coordinates  in  the  cross  section  of  the  left-hand 
waveguide,  and  x'  and  y'  are  rectangular  coordinates  in  the  cross  section 
of  the  right-hand  waveguide.  The  aperture  is  symmetric  about  the  lines 
x  =  a/2  and  y  =  b/2  so  that  the  t^,  t^>  anci  H  directions  introduced  in 
Section  III  can  be  taken  as  the  x,y,  and  z  directions,  respectively. 

The  left-hand  waveguide  is  filled  with  the  loss-free  homogeneous  medium 
(u  e  )  and  the  right-hand  waveguide  is  filled  with  the  loss-free 


Fig.  5.  Two  rectangular  waveguides  coupled  by  a  small 
aperture  in  a  transverse  wall  between  them. 
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homogeneous  medium  (u^,  e^).  It  is  assumed  that  a  >  b  and  that  a'  >  b' 
so  that  the  TE^q  mode  is  the  dominant  mode  in  each  waveguide.  Further¬ 
more,  it  is  assumed  that  the  frequency  is  such  that  only  the  TE^q  mode 
propagates  in  each  waveguide. 

The  excitation  of  the  coupled  waveguides  is  a  TE.^  wave  incident 
in  the  left-hand  waveguide.  If  the  transverse  wall  were  complete  (no 
aperture),  then  in  the  right-hand  waveguide  there  would  be  no  field 
and  in  the  left-hand  waveguide  the  electric  and  magnetic  fields  would 
be  E*a  and  H*a  given  by 


_ia  _  .  -igz  j8z.  .  ttx 
E  =  u  E ,  (e  J  -e  )  sin  — 
—  —y  i  a 


H1*  -  -  u  YE,  («-36z  +  sJ6z)  sin  —  + 
-  —X  i  a 


IttY  ,  -jBz  j8zx  TTX 

u  —  E  (e  -  eJ  )  cos  — 
-z  8a  i  a 


(78) 


(79) 


Here,  u  ,  u  ,  and  u  are  the  unit  vectors  in  the  x,y,  and  z  directions, 
—x  — y  — z 

respectively.  Also,  E^  is  the  amplitude  of  the  incident  wave,  and 


8  = 


Y  = 


1 


k 

a 


n 


a 


where 


k  =  u)  v'v  e 

a  a  a 


n 


a 


(80) 

(81) 


(82) 

(83> 


The  superscript  "a"  in  (78)  and  (79)  indicates  that  the  left-hand 


waveguide  is  taken  to  be  region  a  which  appears  in  the  general  theory. 

The  right-hand  waveguide  will  be  taken  as  region  b  of  the  general  theory. 


In  the  far  zone  (z  <<  0),  the  electric  and  magnetic  fields  £  (M) 

g  A  /v 

and  H  (M)  due  to  M  attached  to  the  left-hand  side  of  the  complete  trans¬ 
verse  wall  can  be  expressed  as 


Ea(M)  = 


„  1  $Z  .  TTX 

u  E  eJ  sm  — 

— y  r  a 


(84) 


Ha(M)  =  -  u  YE  e^Z  sin  — 
—  —  — x  r  a 


jprY  jBz  ttx 
u  -r —  E  e  cos  — 
— z  Ba  r  a 


(85) 


where  E^_  is  an  unknown  coefficient.  Application  of  reciprocity  to  the 
fields  (84)  -  (85)  and  (78)  -  (79)  gives 


Ea(M) 


J  ds  = 


,ia 


(vi-i  + 


v2»2> 


+  E 


ia 


V3M3 


ab 


(86) 


where  S  ^ is  the  cross  section  of  the  left-hand  waveguide  and  J  is  an 

electric  current  sheet  which,  when  radiating  in  the  left-hand  waveguide 

closed  with  the  complete  transverse  wall,  pr  duces  the  electric  and 
is  is 

magnetic  fields  Ji  and 

The  interested  reader  can  show  that  the  current  sheet 


J  =  -  u  2YE,  ejed  sin  —  (87) 

—  — y  i  a 

placed  at  z  =  -d  in  the  cross  section  of  the  left-hand  waveguide  pro- 

is  is  a 

duces  £  and  H  for  -d  <  z  £  0.  In  view  of  the  definitions  of  M^, 

^2 ,  and  M^,  substitution  of  (84)  for  Ea(M),  (87)  for  J^,  (78)  for 
is 

.  .  j  /  n  r\  \  r  .  .  /n  /  \  __ 


33 


E  =  — r 

r  ab 


The  power  P  radiated  by  M  attached  to  the  left-hand  side  of  the 
complete  transverse  wall  is  given  by 


Pa  =  -  Re  Ea(M)  x  Ha*(M)  •  u  ds 


In  view  of  (88),  substitution  of  (84)  and  (85)  into  (89)  gives 


"a  ■  I  I’ll2 


A  development  similar  to  (84)  -  (88)  shows  that  in  the  far  zone 
(z  >>  0),  the  electric  and  magnetic  fields  E^(M)  and  H^(M)  due  to 
M  attached  to  the  right-hand  side  of  the  complete  transverse  wall  are 
given  by 

Eb(M)  =  u  E  e  ^  Z  sin  ^7- 
—  —  — y  t  a 


where 


-  u 

Y’E 

-j(3'z  . 

e  sin 

TTX' 

t 

“X 

t 

a 

u 

jTTY’ 

t  1 

v  0-J3'z 

E  e 

cos 

—2 

B'a 

t 

sf  = 

y? 

-  ^>2 

k.  =  '»>  v'U,  £, 
b  b  b 

n.=/5 


Et  ■  rif 


A  development  similar  to  (89)  -  (90)  shows  that  the  power  P  radiated 
by  M  attached  to  the  right-hand  side  of  the  complete  transverse  wall 


is  given  by 


Comparison  of  (90)  with  (42)  shows  that  the  only  non-zero  element 


of  Y3  is  given  by 


-a  2Y 
hi  ab 


According  to  (98),  the  only  non-zero  element  of  Y^  is  given  by 

vb  =  2 
11  a’b' 


(100) 


It  is  evident  from  (78),  (79),  and  (18)  that  the  only  non-zero  element  of 


the  excitation  vector  is  1^  given  by 


I,  =  2YE . 

1  l 


(101) 


In  view  of  (47),  (48),  (99),  (100),  and  (101),  the  only  non-zero 
element  of  V  of  (58)  is  given  by 


Y  Y’  .  ua  +  ^b 

ab  a'b'  J  4ua  ,1!  U, 

ml  a  b 


(102) 


In  the  far  zone  (z  •'<  0),  the  electric  and  magnetic  fields  due  to 
the  aperture  are  the  fields  due  to  -M  attached  to  the  left-hand  side  of 
the  complete  transverse  wall.  Thanks  to  (84)  and  (85),  these  fields  are 


given  by 


T7a/  w\  17  ™ 

E  (-M)  =  u  E  e  sin  — 

—  —  — v  r  a 


„a,  vr  IbZ  TTX 

H  (-M)  =  u  YE  e  sin  —  + 
—  —  — x  r  a 


(103) 


j"Y  jpz  ”x 
u  “z —  E  e  cos  — 
— z  oa  r  a 


(104) 
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where,  from  (88)  and  (102), 


E  = 
r 


2YE . 

l 


Y  Y ' 

ab  (Ib  +  IMC 


.  l'a  +  ub  , 

J  4“a„lVb 


(105) 


In  the  far  zone  (z  >>  0),  the  electric  and  magnetic  fields  due  to 
the  aperture  are  the  fields  due  to  M  attached  to  the  right-hand  side  of 
the  complete  transverse  wall.  As  given  by  (91)  and  (92),  these  fields  are 


p  ~jB  z  .  -x' 

E  (M)  =  u  E  e  sm  — j- 

—  —  -y  t  a' 


(106) 


Hb(M)  =  -  u 

—  —  -^X 

Y'E 

-j6  Z  .  TX  , 

e  sm  —7-  + 

a' 

u 

— z 

OIL 

F  1*1 

B'a’ 

t  e  cos  r~ 

t  a 

(107) 


where,  from  (97)  and  (102), 


2YE . 


a '  b '  + 


_Y_ 

^ab 


Y' 

a  'b  '  "  J 


U  +  U, 


(108) 


IX.  APERTURE  IN  A  LATERAL  WALL  OF  A  WAVEGUIDE 

In  this  section,  the  moment  solution  is  applied  to  a  rectangular 
waveguide  of  dimensions  a  by  b  coupled  to  a  half-space  by  means  of  a  small 
aperture  in  one  of  its  walls.  As  shown  in  Fig.  6,  the  aperture  is  located 
at  (z  =  0,  x  =  x')  in  the  upper  (y  =  b)  wall  of  the  waveguide.  Except  for 
the  aperture,  the  whole  y  =  b  plane  is  perfectly  conducting.  The  aperture 
is  symmetric  about  the  lines  z  =  0  and  x  =  x'  so  that  the  t^ ,  and  n 


directions  introduced  in  Section  III  can  be  taken  as  the  z,  x,  and  y 


small  f 
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directions,  respectively.  The  waveguide  is  filled  with  the  loss-free 

homogeneous  medium  (y  ,  e  )  and  the  half-space  is  filled  with  the  loss- 

a  a 

free  homogeneous  medium  (y,  ,  e,  ) .  It  is  assumed  that  a  >  b  so  that  the 

b  b 

TE^q  mode  is  the  dominant  mode  in  the  waveguide.  Furthermore,  it  is 
assumed  that  the  frequency  is  such  that  only  the  TE^  mode  propagates  in 
the  waveguide. 

The  excitation  is  a  TE^  wave  which  travels  from  left  to  right  in 

the  waveguide.  If  the  aperture  were  closed  with  a  perfect  conductor, 

then  in  the  half-space  there  would  be  no  field  and  in  the  waveguide  the 

ici  is 

electric  and  magnetic  fields  would  be  E  and  H  given  by 


u 

y 


e 


-jBz 


(109) 


v_  -jBz  .  ttx  JttY  -jBz  ttx 

u  YE.  e  sin  —  +  u  „ —  E,  e  cos  — 
-xx  a  — z  8a  i  a 


(110) 


Here,  E^  is  the  wave  amplitude,  and  B  and  Y  are  given  by  (80)  and  (81) . 
The  waveguide  region  is  taken  to  be  region  a  of  the  general  theory .  The 

half-space  region  will  be  taken  as  region  b  of  the  general  theory. 

a  /* 

Far  from  the  aperture,  the  electric  and  magnetic  fields  E  (-M)  and 
Ha(-M)  due  to  -M  attached  to  the  inside  of  the  upper  wall  of  the  complete 
waveguide  at  (z  =  0,  x  =  x')  can  be  expressed  as 


_a.  r.  „jBz  ,  ttx 

E  (-M)  =  u  E  e  sin  — 

—  —  — y  r  a 


„a.  VT,  jBz  .  ttx  j^Y  jBz  ttx 

H  (-M)  *  u  YE  e  sin  —  +  u  “r —  E  e  r. os  — 
—  —  — x  r  a— zBar  a 


z«0  (111) 


_a.  c.\  „  -IBz  ,  ttx 

E  (-M)  =  u  E  e  J  sin  — 
—  —  -y  t  a 


z»0  (112) 


„a/  r.\  „„  -jBz  .  tx  ,  1_tY  -jBz  ttx 

H  (-M)  =  -  u  YE  e  sin  —  +  u  ^ —  E  e  cos  — 
-  -  - x  t  a  -z  Ba  t  a 


I 
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The  complete  waveguide  is  the  waveguide  with  the  aperture  closed  with  a 
perfect  conductor.  In  (111)  and  (112),  and  E  are  unknown  coefficients. 
Application  of  reciprocity  to  the  fields  (111)  and  (109)  -  (110)  gives 


JJ 

5ab 


Ea(-M)  •  J  ds  =  Hla 


(V^  +  V2M2)  -  E 


ia 


V  M 
3-3 


(113) 


where  S  ,  is  the  cross  section  of  the  waveguide  and  J  is  the  electric 
ab  — 

current  which,  when  radiating  inside  the  complete  waveguide,  produces  the 
electric  and  magnetic  fields  E_  and  H 


It  can  be  verified  that  the  current  sheet 


J  =  -  u  2YE.  e^^  sin  — —  (114) 

—  — y  i  a 

placed  at  z  =  -d  in  the  cross  section  of  the  complete  waveguide  produces 
Eia  and  Hia  for  z>-d.  In  view  of  the  definitions  of  M^,  M,,,  and  , 
substitution  of  (111)  for  Ea(-H),  (114)  for  J_,  (109)  for  E^a,  and  (110) 
for  Hia  in  (113)  gives 


E 

r 


=  J_  rZlUC  v 
ab  0a  1 


+  SV„ 


2 

Ik  S 
a 

6 


V 


(115) 


where 


S  =  sin 


TTx' 


C  =  cos 


7TX 


A  development  similar  to  (113)  -  (115)  shows  that 

i  -i-c  jkis 

Et  ■  lb  '■ sr  V1  -  sv2  -  —  V 


(116) 

(117) 


(118) 


The  power  P  radiated  by  M  attached  to  inside  of  the  upper  wall 
of  the  complete  waveguide  is,  as  calculated  from  (111)  and  (112), 


J 


a  _  abY 
2 


lEr|2  +  i E 1 1  2 ] 


(119) 


Substitution  of  (115)  and  (118)  into  (119)  gives 


a 

P *  V  Y  V 


(120) 


where 


Va  V 

ab 


2  2 
TT  C 

„2  2 

8  a 


0 

TTk2  SC 
a 


TTk2  SC 
a 

e2a 


0 

k4  s2 

a 


(121) 


e2a 


The  half-space  matrix  Y  is  given  by  (63) . 

In  view  of  (47),  (48),  (121),  and  (63),  the  moment  matrix  [Ya  +  Ya] 
reduces  to 

r 

•  0  c 


[Ya+Yb]  = 


0 


(122) 


where 


2  2 
TT  C 


-2 


'1  ,  r,  3,  3Trn,  k  n  a  , 

k  n  3a  b  b  a  a  ml 

a  a 


Vv 
( 2^b 


n  C 2  K  A  u +y. 

,  =  8  S  _ b _ j  f_a _ fcu 

'2  k  n  ab  3Trr,  k  n  m  «,  '  2y, 
a  a  b  a  a  m2  b 


3  2 

k  S 

c«  =  -Arr  + 


^  <^> 


3  n  6ab  3Trn,  r\  a  2c. 


(123a) 


(123b) 


(123c) 


b  e 


C4  = 


~k  SC 
a 

,  2 

n  8a  b 
a 


(123d) 


From  (122) ,  we  obtain 


[Ya+Yb]_1 


-c. 


c  c  —  c 
13  4 


C1C3  "  C4 


-c , 


c  c  —  c 
13  4 


C1C3  “  C4 


In  view  of  (109)  and  (110),  the  excitation  vector  I*  whose  elements 
are  given  by  (18)  reduces  to 


f1* 


(124) 


(125) 


Thanks  to  (123),  substitution  of  (124)  and  (125)  into  (58)  gives 

E  +  E, 


2,,kb  CEiV 


k  aA 
a 


2k? 

t2 


V2  = 


2j3EiSam2 


u_  +  K  qq2  K  Ma 

2(-V±)  +  2J  + 


(126a) 


(126b) 


2n,  E.Sa  M  +  w.  2jkTk  n 

V  - _ b  i  e-  [2  (-2 _ b)  +  >-■*-*  a  ] 

3  n  A  1  V  2y.  ;  3wn.  mlJ 

a  b  b 


(126c) 


where 
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2  2  3 

k.  u  +  u,  c  +  e,  c  +  £,  i  C  k,  k_  rj 

4  ■ 4  r  (~i^— > +  4)<-2f— )(;— r +  3^->  \i 

a  b  b  b  k  3a  b  b 

a 

3  ,3. 


4j 


A  +  ^  ,  ,  (kbkK  ,  kJs2  ,  "V,  (127) 

(  - )  (ft-r—  +  r-r-)  a  +  4a  a  .  tj-)  h— '  +  1“ «rr  +  - H 


2y,  Bab  rj  3iTk  '  e 
b  a  a 


e  ml  3tt  3tttl  k,Bab  '  3. 

b  b  k  Ba  b 

a 


With  the  xyz  coordinate  system  in  Fig.  6,  the  definitions  of 


A  A 

M^»  and  and  formulas  like  (65)  -  (67)  yield 


E  (M,)  =  - 


jkb 


-j  k,r 


2trr 


sin  0  .  , 

- eW 


-jk  r 

. .  D 

b  -  Jkb  e 

E  (M„)  =  (u„  sin  <(>  +  u ,  cos  0  cos  d>)  - z - 

—  —2  —u  — cp  2nr 


2  -Jkbr 

,  k  e 

Eb(M^)  »  -  (uq  cos  6  sin  <p  +  cos  0)  - -  e"1^ 


— 0 


2nr 


where 


(128a) 


(128b) 


(128c) 


ip  =  k,  sin  9  (x'cos  $  +  b  sin  <j>) 
b 


(129) 


Here,  0  is  the  angle  measured  from  the  positive  z  axis,  and  is  the 
unit  vector  in  the  9  direction.  Also,  <p  is  the  angle  measured  in  the 
xy  plane  from  the  positive  x  axis,  and  is  the  unit  vector  in  the  <J> 

a  in 

direction.  Finally,  r  is  the  distance  from  the  origin.  The  factor  e 
in  (128)  is  due  to  the  fact  that  the  aperture  is  located  not  at  the  origin 
but  at  (x=x',  y=b,  z=0) . 

In  the  waveguide  region,  the  far  field  due  to  the  aperture  is  given 
by  (111)  and  (112)  where  and  Et  are  given  by  (115)  and  (118),  respectively. 
The  V's  in  (115)  and  (118)  are  given  by  (126).  In  the  half-space  region 
(y>b),  the  far  field  due  to  the  aperture  is  E^(M)  given  by  (60)  which  is 
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APPENDIX  A.  CONSERVATION  OF  POWER 

In  Appendix  A,  it  is  shown  that  the  moment  solution  for  the  current 

elements  K£  and  l£_  in  Fig.  4  satisfies  conservation  of  power.  Conservation 

of  power  means  that  the  power  entering  the  aperture  from  region  a  is  equal 

to  the  power  flowing  from  the  aperture  into  region  b.  Applied  to  the  situ- 

ation  in  Fig.  4,  conservation  of  power  states  that  the  power  -P  absorbed 

by  the  current  elements  -K£  and  -12.  on  the  left-hand  side  of  the  complete 

screen  in  Fig.  4  should  be  equal  to  the  power  P^  supplied  by  the  current 

elements  K2,  and  (e,  /e  )  l£  on  the  right-hand  side  of  the  complete  screen  in 
—  b  a  — 

3  b 

Fig.  4.  Note  that  the  powers  P  and  P  used  here  in  Appendix  A  are  different 

3  b 

from  those  used  in  the  main  text.  In  the  main  text,  P  and  P  are  calcu¬ 
lated  with  the  impressed  sources  (Jia,  M*a)  and  (Jib,  Mib)  absent.  In 
3  b 

Appendix  A,  P  and  P  are  calculated  with  the  impressed  sources  present. 

The  combination  of  the  left-hand  current  elements  -K2.  and  -I£  is 

equal  to  -M  where  M  is  given  by  (29)  in  which  and  are  unit  magnetic 

current  elements  and  M.  is  the  electric  current  element  12.  =  -jtoe  in  the  n 

— 3  —  a  — 

3  r 

direction.  Hence,  the  power  P  supplied  by  the  combination  of  -K?  and  -12, 
is  given  by 

Pa  *  Re<V1M1  +  V2M2,  Ha(-M)  +  Hla>  -  Re<V3M3>  E3 * * * * * * * II(-M)  +  Ela>  (A-l) 


Thanks  to  (49)  and  (53),  (A-l)  becomes 


a  ia  ia 

P  =  V  Y  V  +  Re<V1M1  +  V^,  H  >  -  RecV^,  E  > 


(A-2) 


By  virtue  of  (44),  the  symmetry  of  Im[Ya] ,  and  the  impulsive  nature  of 
M^,  M2>  and  (A-2)  reduces  to 


a  a-*  *  ia  *  la  *  i3 

P3  =  Re(V  Y  V  +  V.H.  +  +  jut  V.E.  ) 

II  t.  I  a  i  5 


(A-3) 
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The  combination  of  current  elements  K£  and  (e,  Iz  )  II  attached  to 

—  b  a  — 

the  right-hand  side  of  the  complete  screen  in  Fig.  4  is  M  given  by  (29) 

A  A  A 

in  which  and  are  unit  magnetic  current  elements  and  is  the  elec¬ 
tric  current  element  I£  =  -jtoc,  in  the  n  direction.  In  a  development 

b  — 

similar  to  (A-l)  -  (A-3) ,  we  obtain 


Pb  =  Re(V*YbV 


,*  ib  *  ib  ,  *  ib. 

V1H1  "  V2H2  "  ^webV3E3 


(A-4) 


The  sum  of  (A-3)  and  (A-4)  is 

P3  +  Pb  =  Re(V*[[Ya  +  Yb]V  -  I1])  (A-5) 

where  I*  is  the  column  vector  whose  elements  are  given  by  (18) . 

It  is  evident  from  the  moment  equation  (5)  that  the  right-hand 
side  of  (A-5)  vanishes  so  that 

-Pa  =  Pb  (A-6) 

Therefore,  the  moment  solution  for  the  current  elements  K£  and  Ii_  in 

Fig.  4  satisfies  conservation  of  power.  Conservation  of  power  depends 

upon  accurate  evaluation  of  la  and  Re[Ya  +  Yb] .  It  is  interesting  to  note 

s  b  a  b 

that,  as  long  as  Im[Y  +  Y  ]  is  symmetric,  Im[Y  +  Y  ]  has  no  effect  on  the 
right-hand  side  of  (A-5)  and  therefore  no  influence  as  to  whether  power  is 


conserved . 


APPENDIX  B. _ RECIPROCITY 


In  Appendix  B,  it  is  shown  that  the  moment  solution  (3)  for  the 
magnetic  current  M  satisfies  reciprocity.  Reciprocity  states  that 


(E' 


ab 


.la 


.ab 


Mla)ds  = 


,_ba  Tib 
(F.  •  J 


,ba 


Mlb)ds 


(B-l) 


where,  with  reference  to  the  original  problem  shown  in  Fig.  1,  (Eab,  Hab) 
is  the  field  in  region  a  due  to  (Jlb,  Mib) ,  and  (Eba,  Hba)  is  the  field 
in  region  b  due  to  (Jia,  M^a) .  Otherwise  stated,  (E^ab,  Hab)  would  be 
the  field  in  region  a  in  Fig.  1  if  (J^a,  M'*’3)  were  absent,  and  (_Eba,  Hba) 
would  be  the  field  in  region  b  in  Fig.  1  if  (J_^b,  M^b)  were  absent. 

It  is  evident  from  Figs.  2  and  3  that 

(Eab,  Hab)  =  (Ea(-Mb),  Ha(-Mb)>  (B-2) 

(Eba,  Hba)  =  (Eb(Ma)  ,  Hb(M3))  (B-3) 

where  M3  is  the  part  of  M  due  to  (J_la,  Mia) ,  and  Mb  is  the  part  of  M 

due  to  (Jab,  Mib) .  Otherwise  stated,  M3  is  the  value  of  M  which  would 

result  if  (J*b,  M*b)  were  absent,  and  Mb  is  the  value  of  M  which  would 
icl  1 3 

result  if  CJ  ,  M  )  were  absent.  Substitution  of  (B-2)  and  (B-3)  into 
(B-l)  gives 


,„a.„b.  .ia  .ja/.,b.  .,ia.  , 

(E  (M  )  •  J  -  H  (M  )  •  M  )  ds 


(Eb(Ma)  •  Jib 


Hb(MS) 


Mlb)ds  (B-A) 


Equation  (B-4)  is  equivalent  to  (B-l).  Hence,  the  moment  solution  (3) 
for  the  magnetic  current  M  will  satisfy  reciprocity  if  the  moment  solu- 

cl  b 

tions  (3)  for  M  and  M  satisfy  (B-4). 

g  b 

The  moment  solutions  (3)  for  M  and  M  are  given  by 
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Ma  =  vaMt  +  VaM2  +  VaM3 


Mb  =  VbMx  +  VbM2  +  VbM3 


(B-5) 

(B-6) 


&  3.  3 

where  V^,  V2>  and  V3  are  the  elements  of  the  column  vector  V  which 
satisfies 


[Y3  +  Yb]V3  =  lia 


(B-7) 


and  Vb,  Vb,  and  Vb  are  the  elements  of  the  column  vector  Vb  which 
satisfies 

[Ya  +  Yb]Vb  =  Ilb 


(B-8) 


In  (B-7)  and  (B-8),  I^a  and  Iib  are  column  vectors  whose  elements  are, 
with  regard  to  (18) ,  given  by 


4a 


-  H 


ll  '  -  HZ 


ia 

1 

ia 


(B-9) 


and 


4a 


-  j»ca  E3 


ia 


!jb  = 


Tlb  -  Hlb 
T2  H2 


(B-10) 


rib 


!3  =  j^b  E, 


ib 


a  b 

If  it  is  assumed  that  M  and  M  are  given  by  (B-5)  and  (B-6) 
and  that  (J*3,  M*a)  and  0Jib,  M*b)  are  remote  from  the  aperture,  then 
(B-4)  is  equivalent  to 


,„a  .^b.  _ia  ..a-Ab*  ,.Ia*  , 

(E  (M  )  •  J  -  H  (M  )  •  M  )ds  = 


(Fb(M3)  •  Jib  -  Hb(Ma)  •  Mlb)ds  (B-ll) 
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where 


M3  =  +  VaM2  +  V^3 


(B-12) 


aH  b  a  b  a  b~ 

M  =  VXM  +  V2M2  +  V^M 


(B-13) 


Here,  and  $2  are  the  unit  magnetic  current  elements  K£  =  1  in  the 

_t  and  t_2  directions,  respectively.  As  used  in  (B-ll),  the  M3  which 

appears  in  (B-12)  is  the  electric  current  element  I£  =  -juie^  in  the 

n  direction.  As  used  in  (B-ll),  the  M3  which  appears  in  (B-13)  is  the 

electric  current  element  I?.  =  -jwe  in  the  n  direction. 

a  — 

Applying  reciprocity  first  to  the  sources  ftb  and  (J_^a,  M^a)  to 
the  left  of  the  complete  screen  and  then  to  the  sources  Ma  and  (J^b,  Mib) 
to  the  right  of  the  complete  screen,  we  find  that  (B-ll)  is  equivalent  to 


Vb  H*a  +  V^H2a  +  jweaV3E3a  -  -  VaH^b  -  VaH2b  -  jUi^E^ 


(B-14) 


Equation  (B-14)  is  always  equivalent  to  (B-ll)  regardless  of  the  values 
of  the  V's. 

It  will  now  be  shown  that  the  V's  which  satisfy  (B-7)  and  (B-8) 
also  satisfy  (B-14).  In  matrix  form,  (B-14)  is 

(B-15) 


-vbila  =  -  vaiib 


Va  and  Vb 


-  b 

where  V  and  V  are  the  transposes  of  the  column  vectors 
in  (B-7)  and  (B-8) .  The  column  vectors  I*a  and  l^b  in  (B-15)  are  the 
same  as  those  in  (B-7)  and  (B-8).  Multiplication  of  (B-7)  by  Vb  from 
the  left  gives 

Vb[Ya  +  YblVa  =  Vbfia  (B-16) 

Multiplication  of  (B-8)  by  V3  from  the  left  gives 

Va[Ya  +  Yb]Vb  -  Vatlb 


(B-17) 
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Because  [Ya  +  Y^]  is  a  symmetric  matrix,  the  left-hand  sides  of  (B-16) 

and  (B-17)  are  equal  to  each  other.  Therefore,  the  right-hand  sides 

of  (B-16)  and  (B-17)  are  equal  to  each  other  so  that  (B-15)  is  true. 

In  the  previous  paragraph,  we  found  that  the  column  vectors 

Va  and  of  the  V's  in  the  moment  solutions  (B-5)  and  (B-6)  satisfy 

(B-15).  Backtracking,  we  see  that  these  V's  satisfy  (B-14)  which  is 

equivalent  to  (B-ll) .  Since  (B-ll)  is  equivalent  to  (B-4) ,  we  con- 

3  b 

elude  that  the  moment  solutions  for  M  and  M  satisfy  the  statement 
(B-4)  of  reciprocity.  In  other  words,  the  moment  solution  (3)  for  the 
magnetic  current  M  satisfies  reciprocity. 
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